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£> . Abstract 

o 

The goal of this article is to give a formal derivation of Ohm's law of Magnetohydro- 
dynamics (MHD) starting from the Vlasov-Maxwell-Boltzmann system. The derivation 
is based on various physical scalings and the moment methods when the Knudsen num- 
ber goes to zero. We also give a derivation of the so-called Hall effect as well as other 
limit models such as the Navier-Stokes-Maxwell system. Our results include both the 
Oh ' compressible and incompressible MHD models. 
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1 Introduction 

There are different models to describe the state of a plasma depending on several parameters 
such as the Debye length, the plasma frequency, the collision frequencies between the different 
species. One of the goals of this project is give a relation between the several models used, 
namely how to derive rigorously the simpler models from the more complete ones and find 
the regimes where these approximations are valid. Formal derivation of these models can be 
found in Plasma Physics text books (see for instance Bellan [9], Boyd and Sanderson |12| . 
Dendy p2] and the paper [10] etc.) 

Indeed, since the plasma consists of a very large number of interacting particles, it is 
appropriate to adopt a statistical approach. In the kinetic description, it is only necessary 
to evolve the distribution function f a (t,x,v) for each species in the system. Vlasov equation 
is used in this case with the Lorentz force term and collision terms. It is coupled with the 
Maxwell equations for the electromagnetic fields. 

If collisions are important, then each species is in a local equilibrium and the plasma is 
treated as a fluid. More precisely it is treated as a mixture of two or more interacting fluids. 
This is the two-fluid model or the so-called Euler-Maxwell system. 

Another level of approximation consists in treating the plasma as a single fluid by using the 
fact that the mass of the electrons is much smaller than the mass of the ions or by using that 
collisions yield that all particles will evolve (at leading order) with the same macroscopic 
velocity. This yields the Hall-MHD (Hall magneto- hydrodynamic model). Then, one can 
derive MHD models if the Hall effect is negligible. 
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The question on how the kinetic theory and the fluid dynamics are related is very inter- 
esting both physically and mathematically and it goes back to the founding work of Maxwell 
[40] and Boltzmann [11]. Moreover, the purpose of the Hilbert's sixth problem [30 1 is to seek 
a unified theory of the gas dynamics including various levels of descriptions from a math- 
ematical point of view. There has been a lot of important progress on the hydrodynamic 
limits from the Boltzmann equation over the years; so far there are essentially three dif- 
ferent mathematical approaches. The first is based on spectral analysis of the semi-group 
generated by the linearized Boltzmann equation; see [7] [33] [41]. The second is based on 
Hilbert or Chapman-Enskog expansions [13[ [16], [24"] ; see more recent work in [26[ [28] [29] [32] . 
The third approach, initiated by Bardos-Golse-Levermore [H [5], is working in the frame- 
work of DiPerna-Lions' renormalized solutions [18] , to justify global weak solutions of incom- 
pressible flows (Navier-Stokes, Stokes, and Euler), and (compressible) acoustic system; see 
[3[g[l9l[22l[3ll[35l[36l[38l[l2] (see also [2] for the non-cut-off case). 

However, there are only limited results on hydrodynamic limits for charged particles, 
which satisfy the Vlasov-Maxwell-Boltzmann (VMB) system, even at the formal level due 
to the complexity of the system and its underlying muti-scale feature. One rigorous result 
is given in [31] where a diffusive expansion to VMB system was studied with one particular 
scaling in the incompressible regime in the framework of classical solutions [25J, and as a 
by-product, new fluid equations interacting with the electric field, where the magnetic effect 
appears only at a higher order, were derived. See also [8] [27] [39] for hydrodynamic limits 
from the Boltzmann equations in the presence of an interacting field, but without a magnetic 
field. 

In this article, we are interested in the derivation of MHD type equations, which describe 
the motion of electrically conducting media in the presence of a magnetic field, from the VMB 
system by introducing various scalings. As a result, we can also identify the corresponding 
Ohm's law of the MHD equations. The goal here is to provide formal expansions that show 
the relevant scalings. Our analysis is the first step towards a more rigorous analysis. It can 
be seen as the extension of the paper of Bardos-Golse-Levermore [4] to the VMB case. Some 
rigorous derivation will be given in a forthcoming paper. 

Let us also mention that various MHD models were obtained from the (macroscopic) 
two-fluid Euler-Maxwell equations in [10] by taking different asymptotic limits. Here, our 
goal is to start from the kinetic level. 

1.1 Two species Vlasov-Maxwell-Boltzmann system 

The dynamics of charged dilute particles (e.g., electrons and ions) is described by the Vlasov- 
Maxwell-Boltzmann system: 

d t F + + v ■ X7 X F + + ^±-(E + v x B) ■ V V F + 
m + 

d t F- + v ■ V X F- - — (E + v x B) ■ V„F_ 

m_ 

F±(0,x,v) = F 0> ±(x,v). 

Here F±(t,x,v) > are the spatially periodic number density functions for the ions (+) and 
electrons (-) respectively, at time t > 0, position x = (xi,X2,x^) G [ — 7r,7r] 3 = T 3 , velocity 
v = (vi,V2,Vs) £ R 3 , and e±, m± the magnitude of their charges and masses. 

The collision between particles is given by the standard Boltzmann collision operator 
Q(G±, G2): Let G\(v), G2(v) be two number density functions for two types of particles with 



= Q(F+,F + ) + Q(F+,F_), 

= Q(F_,F+) + Q(F_,F_), (1.1) 
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masses mi and diameters <7j (i = 1,2), then (p 83 and p 89 in [TJ]) Q(Gi,G2)(v) is denned 
as 

7 On +cj 2 ) 2 [ b{u-v,uj){G 1 {v)G2{u)-G 1 {v)G2{u)}dudu (1.2) 
4 7R3x5 2 

= <5gain(Gl, G2) — Q\oss{G\, G2). 

For hard-sphere interaction, the collision kernel b is given by b{u — v, uS) = \{u — v) ■ w\ (see 
|14j ) and in this article, we assume the hard-sphere interaction, but the formal derivation 
will be valid for other general collision kernels. Here weS 2 , and 

' 2m2 \( \ 1 / 1 2mi 

f = t> \(V — U) ■ U!\UJ , u = u -\ \(v — u • UJ \UJ. 1-3) 

mi+m 2 mi+m 2 v ; 

which denote velocities after a collision of particles having velocities t> , u before the collision 
and vice versa. Notice that the elastic collision (|1.3|) implies the conservation of momentum 
miv+m,2U and energy ^mi|v| 2 -|- ^m,2|u| 2 during the collision process. To clarify the collisions 
between two types of particles, we use the following notation: 

Q+ ^Q(F + ,F + ); Q± = Q(F + ,F_); ^ Q(F_, F + ) ; Q~ ^ Q(F_, F_) . 

Note that Q + and Q~ are the usual collision operators of one species. 

The self-consistent, spatially periodic electromagnetic field [E(t, x), B(t, x)] in (11. ip is 
coupled with F(t, x, v) through the Maxwell system: 

fio^odtE — V x B = —fioJ = —^0 / V i e +E+ ~ e-F_}dv, 

d t B + V x E = 0, V • B = 0, 

V • £ = — / {e + F + - e-F_}dv, 
eo Jr3 

£(0,x) =£ (x), £(0,z) =BoOzO- 

Here eo and are called the electric permittivity and the magnetic permeability of the 
plasma. And the speed of light c is given by 



It is well-known that for classical solutions to the VMB system, the following conservation 
laws of mass, total momentum (both kinetic and electromagnetic) and total energy (both 
kinetic and electromagnetic) hold: 



d_ 

Tt 

d_ 

dt 7ts x r3 



/ m + F + (t) = 0, 

JT 3 xR3 

m_F_(t) = 0, 

yAI v(m + F + (t) + m-F-(t)) + ]- [ E(t) x B(t)\ = 0, 

"* UT 3 xR 3 47T J T 3 J 

\\! H 2 (m + F + (i) + m_F_(i)) + ^- / |tf(t)| a + |B(t)| 2 )=0. 

l2J T 3 xR 3 87T J T 3 J 



_d f 1 
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Moreover, we also have the celebrated H-Theorem of Boltzmann 

i{/ (F + (t)lnF + (t) + F„(t)lnF„(t))\<0. (1.4) 
dt Ut 3 xr 3 J 

Indeed, the above conservation laws and the entropy dissipation can be derived from 
the following well-known properties of the collision operators Q + , Q^, Q^, Q~ . We refer to 
[H EJ IS] for more details. 

(i) Mass conservation: 

J Q+dv = , j Q ± dv = , j Q T dv = , J Q~dv = . 

(ii) Momentum conservation: 
Ctm+vdu = , f(Q±m + v + Q^v)dv = , / QTm-vdv = . 

(hi) Energy conservation: 

/ <^ M '* = , /(«^M' + Q*m-M')* = , / fl-m- W* = ■ 

(iv) Entropy inequalities: 

J Q + {F + ,F + ) In F + dv < 0,J Q~(F-,F-)]nF-dv <0, 

J (Q ± (i ? +,i ? -) In F+ + Q T (F_,F+) In FJ)dv < 0. 

As a consequence of the entropy inequalities, we can find distribution functions (two 
Maxwellians) which cancel the collision operators. 

(v) Local Thermodynamical Equilibria: 

M + (t,x,v) = n + (^) 3/2 e- m +^"l 2 / 2T , M.(t,x,v) = n _^f/* e -™-\^ /*T (L5) 
so that 

Q + {M + ,M + ) + Q ± (M+,M„) = and Q T (M_ , M + ) + Q~ (M_ , M_ ) = 

where n+ and n_ are the ion and electron density and u and T are the common mean velocity 
and temperature and they may depend on t and x. When the macroscopic variables do not 
depend on t and x, we call them global Maxwellians. When u = 0, we define 

M«) = n + (^) 3 / 2 e-+H 2 /2T 0) M » = n_(^) 3/2 e- m -H 2 / 2 T . (L6) 

Note that when e+n+ = e_n_, the global Maxwellians and H-{v) together with 

E = B = define a stationary solutions to the VMB system, while the local Maxwellians are 
not necessarily solutions. 

We now introduce the linearized collision operators over the velocity space for Q + + 
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and QT + Q around M+ and M_ in the following vector form: 

.A 

M^{Q+(M + f + , M+) + Q+(M + ,M + f + ) + Q ± (M + f + , M_) + Q ± (M + , M_/_)}\ (L7) 
MZ l {Q T (M_ /_ , M + ) + (M_ , M+/+) + Q- (Af_ /_ , M_) + Q~ (M_ ,M./_)}j' 

We define the following inner product in t>: 

(/:)'C:))M^/ R . (M+ ^ ++M - / - 9 - ,<fc 

and denote by L 2 M the associated Hilbert space. We will use (-,•) to denote the usual L 2 
inner product without the weight. We summarize the properties of £ in the following lemma. 
The proof can be found in [I] (pp 635-638). 

Lemma 1.1. We assume the hard-sphere interaction for the collision kernel. 
1. £ is the sum of a diagonal operator f — > vf 

vf= ( U+( i"\ f f + ) with i/±(H)~l+H 
and a compact operator A. The domain of £ is given by 

£>(£) = {/: ||(1 + H)I/|U <oo}. 



2. £ is self-adjoint in L\,: 



M ■ 

f+\ /(f+\ n (g^ 



f-J 1 \q-J/m \\f-)' V .•'/- J i i ' 

£ is non-negative. 

3. The kernel of £ is a six- dimensional linear space: 



ker £ = Span 



2 

m-Vij \m—\v\ 



0\ fm+Vi\ { m+\v\ 
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4- Any function f £ D(£) can be written as f = qj+Wf withqf G ker£ andwj S (ker£) 
and we have (<£/,/) m > ^o||(l + |v|) 2 tf/|| 2 - 

Note that this Fredholm operator £ can be inverted after checking that the inhomogeneity 
is perpendicular to its six-dimensional null space 



1\ /0\ (m+v A (m+\v 



0/ VI/ \m-ViJ \m-\v 



2 
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{4> , ...,(p 5 }- 



1.2 F-G formulation 



In the MHD model, the medium is considered as a single fluid, in other words, one neglects 
the difference in motion of the electrons, various kinds of ions, and neutral particles. To 
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derive the MHD equations from the VMB system (jl.ip . it is convenient to consider the total 
mass density and the total charge density for F + and F~: 

F = m + F+ + m_F_ and G = e+F + - e_F_ . (1.8) 

Then the VMB system (jl.ip can be rewritten as 

d t F + v ■ V X F + (E + vxB)- V V G = m + (Q + + Q ± ) + m_(Q T + Q~), 

d t G + v ■ V X G + (E + v x B) ■ V V {-^^F + e + m - ~ e ~ m + G} (L9) 

= e+(Q + + Q ± ) - e_(Q T + Q~), 



and 



H e dtE - V x B = -/x / u Gcfo, 9 ( B + Vx£: = 0, 
V • E = — [ Gdv, V • B = 0. 



(1.10) 



Notice that the Maxwell's equations are coupled with only G. We will use both F + ,F~ 
formulation (|1.1|) and F, G formulation (|1.9|) throughout the article. On one hand, we remark 
that when m+ = m_(= m), the collision operators are of one-species and thus the right-hand- 
sides of (|1.9|) can be significantly simplified: 

• m + (Q + + Q±) + m_(Q T + Q~) 

= m{Q(F + ,F + ) + Q(F+,F_) + Q(F_,F + ) + Q(F_,F_)} 
= mQ(F + + F_ , F + + F_ ) 

= -Q(f,f) 

• e+(Q + + Q ± )-e_(Q T +Q-) 

= e+{Q(F+,F + ) + Q(F+,F_)} - e_{Q(F_, F + ) + Q(F_ , F_)} 
= Q(e+F + -e_F_,F + + F_) 

= -Q(G,F) 
m 

In this case, the linearized operator (II. 7ft is orthogonally decomposed: the linearized operator 
of Q(F, F) has the five-dimensional kernel {1, v, \v\ 2 } and the linearized operator of Q(G, F) 
has the one-dimensional kernel {1}, for instance see [8j [31] for the use of these linearized 
operators. This special setting will be frequently used. On the other hand, we remark that 
when e+m__ = e_m+, then the term e+ "^~^ l ~ m+ G cancels in (|1.9p . 
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2 Near local thermo dynamical equilibria 



By using the collision invariants of the Boltzmann collision operators, if F, G, E, B satisfy 
the VMB system (|1.9p and (jl.lOp . we formally deduce the following local conservation laws: 

• dt( J Fdv^j + V X -(J vFdv} =0 

Fdv) = ( I Gdv)E+ ( I vGdv) x B 

(2.1) 



• dt{ J vFdvj + V 

• dt( I ^-Fdv) +V X ( / r^-Fflr) F- { / cGdc 



2 J V J 2 

• fl^ y Gcfo) + V x -( K j vGdv} =0 

which describe the conservation of the total mass, momentum, and energy and the conserva- 
tion of the charge density. Of course, the equations are coupled with (|1.1U|) . However, the 

integrals inside Va;- such as J v ® v Fdv, fv^- Fdv J v Gdv are not functions of the other 
macroscopic quantities and they depend on the entire distributions F and G unless the ap- 
propriate forms for F and G are assumed. Thus, in general, the moment method does not 
give a closed form of a finite set of equations and this is often referred as the closure problem. 
In order to close the system of balance laws (|2.ip , we look for the solutions very close to local 
thermodynamical equilibria F + ~ M + and F^ ~ M_ such that in such asymptotic regime, 
the above integrals can be expressed in terms of the given macroscopic variables. This ap- 
proximation can be realized by rescaling the system, namely by comparing the size of the 
system and the time scale of interest with the typical space and time scales of collisions. The 
task is not all trivial because the physical regimes are much richer than the pure Boltzmann 
equation and one needs to deal with multi-scales encoded in the VMB system. 

In what follows, we introduce a few different hyperbolic scalings to the VMB system 
and study their hydrodynamic limits formally. Also, we discuss the Hall effect of the MHD 
equations. We want to point out that we did not try to get all possible scalings. In particular, 
one can easily derive simplified models to those we have here. 

2.1 l|— fluid Euler-Maxwell system with the common velocity and tem- 
perature 

We start with the simplest possible hyperbolic scaling to (II. ip : Let t = st, x = ex, F± = ^F±, 
b = e 2 b, E = -E and B = \B. Dropping ~, we obtain the following rescaled VMB system: 

d t F% + v ■ V x Fl + ^(E e + vxB £ )- V v Fl = ±{Q+(F*,F%) + Q ± {F" + , Ft )}, 
m_|_ £ 

d t F e _ + v ■ V x Fi - — (E £ + v x B £ ) ■ V V F1 = -{Q^(Fi,F £ + ) + Q~(Fi,Fi)}, 

m_ e 

Hoe d t E £ - V x B e = / v{e+F £ + - e-Fl}dv, d t B e + V x E £ = 0, 

V-E £ = —[ {e+Fi - e^Fljdv, V-B £ = 0, 

where e > is a Knudsen number. 



(2.2) 
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Suppose F± — )■ F±, [E £ ,B £ ] — > [E,B] as e — > 0. Then from the first two equations in 
(12. 2|) . we see that 



Q+ ( F + ,F + ) + Q ± (F + ,F-) = 0; QT(F-,F + ) + Qr [F~ , F~ 
from which we deduce that 







F 



M+ and F~ 



Al- 



locs! Maxwellians given in (|1.5p . There are six fluid variables n + , n„, u and T to be deter- 
mined: the dynamics of n+, n_, u and T can be derived by taking the moments of the first 
two equations. The system of local conservation laws 



d t F% + v ■ V X F £ + + ^±-{E £ + v x B £ ) ■ V V F: 
8 t Fl + v ■ V x Ft - ^-(E £ + v x B £ ) ■ V V F_ 



0, i = 0,...,5 



(2.3) 



is, in general, not closed. If each moment converges, we can pass to the limit of e — > 0. To 
write the equations for n+, n_, u and T, we first compute the first few moments of local 
Maxwellians M±: 



M±dv = n±, 



vM±dv = n±u, 



\v\ 2 M±dv = n±\u\ 2 + 3^±^, 

m± 



v ® vM±dv = n±u (8) u + 
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IfpvM-tdu = n-t|u| 2 u + 5 



n±T 
m± 



-u. 



Thus from (12.3|) where i*± are placed by i<± 



<9 t n+ + V x -(n + u) = 
<9tn_ + V x ■ (n_u) = 

<9f (m + n + + m_n_)u 



(m + n + + m_n_)u ® u 



+ V, 



(n + + n_)T) 



= {n + e + — ri-.eJ)E + (n + e + — n_e_)u x 
d t 



(2.4) 



1 3 

-{m + n + + m_n_)|u| 2 + ^Oh + n -)T 



1 5 

-(m_|_n + + m_ra_)|u| 2 u + + n_)Tu 



(n + e + — n_e_)u 



with 



HoeodtE — V x B = — /xo(e+n+ — e_n__)u; c\-B + VxE = 0; 

V • E = — (e+n+ - e_n_); V • B = 0. 

By introducing the total density p and the charge density a: 

m + n + + m_rt_ =: p and e + n + — e_n_ =: a 



(2.5) 
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the above system (|2.4p and (|2,5p can be equivalently written as 



• d t p + V x • (pu) = 

• d t (pu) + V x -(pu®u) + V x (ra + + n_)T = 




(2.6) 



+ -(n + + n_)Tu = aE-u 



d t a + V x -{au) =0 



where 



p eod t E -V x B 



poau ; 9 t B|Vx£ = 0; V ■ E = —a ; V ■ B = (2.7) 



and 



and n + + n 



(e + + e_)p + (m_ — m + )a 



(2.8) 



ra± = 



e_m + + e + m 



e_m + + e + m_ 



The first three equations in (I2.6P represent the balance laws of total mass, total momentum, 
and total energy and the last equation is the conservation of the charge density. In this 
regime, the current density J is given by au. Hence, by letting 



the solvability condition (|2.3[) yields Euler-Maxwell type equations (|2.6p and (|2.7p having the 
same mean velocity and temperature fields which are therefore not exactly of two-fluid. From 
(|2.8p . note that for a special case of m + = ni-, the pressure is induced by the total density 
but in general, it depends on the charge density as well. 
We have proved the following theorem. 

Theorem 2.1. Let F±, [E e ,B £ ] be the classical solutions to (12. 2p . Assume that F± — > F±, 
[E £ ,B £ ] — > [E,B] strongly as e — > 0. Then the limit F± is a Maxwellian given by (jl.5p and 
moreover, n±,u,T, E,B solve the l^-fluid Euler-Maxwell system (12. 4p and (12. 5ft . 

For a rigorous proof of this result, we can argue as in Caflisch [13]. Indeed, the derivation 
given above can also be used to provide a Hilbert expansion as in [13J. This is then used to 
provide an error estimate between the solution to the 1^-fluid Euler-Maxwell system and the 
solutions to the VMB system (I2.2p in the limit when e goes to zero. 

2.2 Compressible resistive MHD 

We have seen that it is convenient to write the equations in terms of the total density and the 
charge density in the previous section. In this subsection, we start with the F—G formulation 
introduced in Section fl.21 We set m+ = m_ = e+ = e_ = 1 for simplicity of the presentation. 
Then by using (jl.llj) . the VMB system (II. 9h reads as follows: 



In order to see the ideal MHD type equations with Ohm's law, we consider the quasi-neutral 
regime so that G ~ 0. We next introduce another scaling to (|2.9p as follows: Let t = et, 



x = ex, E = XE, B = 4.B, G = -^G, F = ~F, and b = eb where E = 0(y/e), 



p = (n + + n_)T (ideal gas law), 



d t F + v ■ V X F + (E + vxB)- V V G = Q(F, F), 
8 t G + v V X G + (E + v x B) ■ V V F = Q(G, F). 



(2.9) 
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B = 0{y/e), G = 0(e 3 / 2 ) and F = O(e). Dropping ' ", we obtain the following rescaled VMB 
system: 



d t F e + v ■ V X F £ + {E £ + v x B £ ) ■ V V G £ = ^Q{F £ , F £ ), 
d t G e + v ■ V X G £ + + V -^-) ■ V V F £ = - £ Q{G e , F £ ), 
coupled with the Maxwell's equations 

Ho€ d t E £ - VxB E = -/i J vG £ dv, d t B e |VxF = 0, 

V • E £ = — [ G £ dv, V • B e = 0. 



(2.10) 



(2.11) 



Suppose that there exists a classical solution [F £ ,G £ , E £ , B £ ) to (|2.10p and (|2.1ip and that 
there exist F, G, E, B so that F £ ->■ F, G £ -)• G, E £ -)• E, B £ -> 5 as e -> 0. Then from 
(|2.10p , as e — )• we deduce that 

Q(F,F) = F = M = — Pe-W/xr ( 2 .12) 

also 

(£ + wxB)-V/ = Q(G,F). (2.13) 

This case, it is already different from the previous ones in that G is not determined by a 
Maxwellian. Of course here G is a higher order fluctuation. Define 

Cg = -^Q(Mg,M) 

where M is the Maxwellian given in (|2.12p . Then from Lemma 1 1 . X I we deduce that ker£ = 
span{l}. The solvability condition for G is automatically satisfied and thus by inverting the 
linearized operator C in (|2.13p we derive that 

G = -M[£- 1 {—(E + vxB)-V v M} + iJja], where = 0. (2.14) 
Let us define the first two moments of G by a (the charge density) and J (the current): 

a ■= J Gdv and J := J vGdv. (2.15) 
Then from (|2.14p we can write the current J as 

J = j vGdv = (v,G) = -(^Mr^^ + jix B)-V V M}) 



and since 



j^(E + v x B) ■ V V M = -{E + v x B) ■ =~(£ + uxfl).(»-u) 
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and by using the fact that C 1 acts only in v, we obtain 

J = (v, MC- X {^[E + u x B) ■ (v - u)}) = (vM, ~(E + u x B) ■ C' 1 ^ - u)> 
= ((« - u)M, ^(i? + ux5)- - u)). 

Since J" (?/ — u*)vC _1 (z;- ? — u^)Mdv = for i 7^ j, the last expression for J can be reduced to 
the following: 

J =-(E + ux B) (2.16) 
V 

where 

1 1 f 

- = — (v -u) ■ C^(v - u)Mdv > 0. 
rj 3T y 

Here, the number 3 represents the space dimension 3. The relation (|2.16p is called the Ohm's 
law. rj is called resistivity and 1/77 is the conductivity. Next we take the moments of F 
equation in (|2.10p . then by collision invariants of Q, we get 

{dfF e + „ . VxF e + (E £ +VX B e ) ■ V v G e , {l,v, 1^}) = (2.17) 



and 



fi £od t E £ - V x B £ = -no I vG £ dv, d t B £ + V x E £ = 

V • E £ = — I G £ dv, V • B £ = 0. 



(2.18) 



We pass to the limit and by using the fact that F is a Maxwellain for which we can explicitly 
compute the moments and by using the notation introduced in (|2.15p , we obtain the following 
resistive MHD type equations for p, u, T, a, J, E, B: 

• d t p + V x -{pu) = 

• dt(pu) + V x -(pu(g)u) + V x {pT) = oE + J x B 



dA P ^ + lpT\ + V x .(^ + lpTu)=E.J 



(2.19) 



2 2' J V 2 2 

• d t cr + V x -J = 

• J= -(E + uxB) 

V 

• p e d t E - V x B = -no J, d t B + V x E = 0, V • B = 0, V ■ E = — . 

eo 

We have proved the following theorem. 

Theorem 2.2. Let [F £ ,G £ ], [E £ ,B £ ] be the classical solutions to ([2TTU]) and (|23T|) . Assume 
that [F £ , G £ ] -4 [F, G], [E £ , B £ ] ->• [E, B] strongly as e ->• 0. T/ien f/ie Zimzi F is a Maxwellian 
given by ()2.12p and moreover, p,u,T, o~,J, E,B solve the electromagneto-hydrodynamic sys- 
tem (I2TT9D . 

The electromagnetic field in the system (|2.19p is still governed by the full Maxwell's 
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equations; if we consider the regime where the electric permittivity eo is sufficiently small, 
namely if we consider the quasi- neutral limit, then we would get a new system in which 
the magnetic field becomes dominant. We may set eo ~ e by fixing po i n the Maxwell's 
equations (I2.1ip to obtain different approximations to the dynamics of the electromagnetic 
field. Namely, if we consider 

Ho£d t E £ - V x B £ = -(iq J vG e dv and V • E e = i J G £ dv, 

following the same procedure, we would get V x B = p^J and a = 0. And (12.190 reduces to 

• d t p + V x -{pu) = 

• d t (pu) + V x -(pu®u) + V x {pT) = J x B 

. at (^ + \ P T) + ^(^ + \ P Tu)-E.J (2 , 0) 

• J= -(E + u x B) 

V 

• V x B = no J, d t B + V x E = 0, V • B = 
which are the well-known compressible resistive MHD equations. 

Theorem 2.3. Let [F £ ,G £ ], [E £ ,B £ ] be the classical solutions to (pHO]) and IpTT]) where 
eo is taken to be e. Assume that [F £ ,G £ ] —?- [F,G], [E £ ,B £ ] [E,B] strongly as e — > 0. 
Then the limit F is a Maxwellian given by (12.120 and moreover, p, u, T, J and E, B solve the 
compressible resistive MHD equations (|2.20p . 

We remark that the ideal MHD equations are valid when the resistivity r/ is negligible. 

Also, we point out that the simplification of setting the physical constants to be one allows 
us to treat the collision operators decoupled as in (|2.9p . The General case, in particular when 
m + ^ m_, can be treated similarly. In such case, we have to start with (jl.ljl or (11.90 and 
need to treat the collision operators as matrix form as in (jl.7p and use that it has a six 
dimensional kernel. We do not detail this here. 

2.3 When e_m + ^ e + m and Hall effect 

The Hall effect is known as two-fluid effect. In our setting, the Hall effect can be captured 
when e_m + ^ e+m_ at the higher order. For instance, we take the same scaling as done in 
(12~T0|) to CCS]), 

d t F £ + v ■ V X F £ + (E £ + vx B £ ) ■ V V G £ = ^{m + (Q + + Q±) + m_(QT + q-)^ 

„ _ „ e ,E £ v x B £ . _ t e+e- „ i e + m_ — e_m + 

d t G £ + v ■ V X G £ + ( — + ) • V„{— F £ + ea -t ±G £ } 

e e m + m- m + m- (2 21) 

> v ' v ' ' 

(*) 

= - £ [e + {Q + + Q ± )-e-{QT + Q-)]. 

The main difference is the presence of the term (*), which will affect the Ohm's law. 

Following the same spirit of the Chapman-Enskog expansion, if we also keep the 0(y/e) 
terms in the second equation of (I2.2ip . instead of (12.16p . we would get the following Ohm's 
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law 

i 

v j = ^—( E + uxB) + £ -l{ e+m -- e ~ m+ ){aE + JxB) (2.22) 

where J x B is often referred to as the Hall effect term. We remark that this Hall effect is not 
seen by the moment method that we have used for the previous theorems but rather this can 
be captured by the higher order Chapman-Enskog type expansions. It would be interesting 
to justify this higher order effect in a rigorous framework. 

Remark 2.4. We have chosen our sailings so far based on the appropriate time, space and 
forcing with fixed masses m + and m_ . When m_ <C m + and is comparable to Knudsen 
number, one can also parametrize ^j—- By doing so, for example, we can get the rescaled 
system (|2.10p or (|2.2ip without rescaling forcing terms. In such regime, one may want to 
study the appropriate collision operators as suggested in [15]. 



3 The incompressible regime 

As we have seen in the previous section, the derivation of hydrodynamic equations from the 
Boltzmann equation is closely related to finding approximate solutions of the Boltzmann 
equation, because the fluid variables are defined and change in space and time scales that 
are very large when measured in the units of the mean free path and mean free time between 
collisions. In the present section, we'd like to discuss the longer time scaling. We will mainly 
work on the simplified case where the constants are taken to be unity and present general 
cases at the end when we discuss the Hall effect. 

3.1 Incompressible Navier-Stokes-Fourier-Maxwell equations 

Recall the system (|2.9|) with the Maxwell's equations: 

8 t F + v ■ V X F + (E + v x B) ■ V V G = Q(F, F), 
8 t G + v V X G + (E + v x B) ■ V V F = Q(G, F), 

HQtQd t E - V x B = -/j, / vGdv, V • B = 0, (3-1) 

d t B + V x E = 0, V • E = — [ Gdv, 

where F = F + + F_ is the density for the whole particles and G = F + — F_ represents the 
disparity between two species (the charge density). 

We will discuss how to obtain the incompressible Navier-Stokes-Fourier-Maxwell system. 
To keep the hyperbolic structure of Maxwell's equations, we take and eo as 0(1) and let 
t = et, x = ex, v = ev (longer time scale), E = eE, B = eB, and F = eF. Then the system 
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(|3.ip can be rescaled as follows: by dropping ~ 

ed t F e + v ■ V X F £ + (eE E + v x B £ ) ■ V V G £ = \q{F £ , F £ ), 

e e z e 
H e dtE £ - V x B £ = -no [ v G £ dv, V ■ B £ = 0, 

<9 t -B £ + V x £ e = 0, V • E £ = — [ G £ dv. 

We now consider the following perturbation of F £ ,G £ ,E £ ,B £ around \p, 0,0, 0]: 

F £ = + eft), G £ = eufi, E £ = eEf, B £ = eB{ 
where fi is a global Maxwellian normalized as 

1 -M72 



(3.2) 



(2vr)3/2 



(3.3) 



Suppose that there exist the classical solutions ff, g\, Ef, B\ to (|3.2|) and further assume 
that as e — > 0, ff — > f, g\ — > g, E\ — > E, B\ — > B. Under this convergence assumption, we 
will derive the dynamics of the limiting variables f,g,E,B. We first rewrite the system 
in terms of the perturbed variables: 

ed t fl + v • V x ff + e(eEl + v x B\ ) • -V v (wf) 

A* 

= j-{Q(l*fl,l*) + Q^fJLf!) +eQ{fifl^fi)} 

E £ v x B £ 1 
ed t g{ + v ■ V x gf + (— + =— ) • -V„(m + e/x/f ) 

= — {Q(«rf, AO + )}, 
HoeodtEl - V x Bf = -/j, / u/i#f dv, V • Sf = 0, 

d t B{ + V x Ff = 0, V • E\ = — [ Hgfdv. 

We recall the linearized collision operators L and C (see for instance [8j [31]) : 

Lf = --{Q(fj,f,fj,) + Q(fj,,fj,f)} and Cg = --Q(fig, fi). 

Note that is equivalent to the linearized collision operator introduced in (jl.7p . We also 

recall that kerL = span{l,i>, \ v | 2 } and ker£ = spanjl}. 

By taking e — > in the first two equations in fj3.3j) and using the convergence assumption, 
we first obtain 

= —Lf (3.4) 
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and 

-E-v+-(yxB)-V v (jif) = -Cg. (3.5) 



Since / is in the kernel of L by (13. 4p . / can be written as 
We also introduce 

a := gpdv and J := gvpdv. 



From (13.51). we can write the current J as 



J = J gvfidv = J C 1 (E ■ v — ~(v x B) ■ V v (fJ,f))v/idv. 

Since 

(v x S) • V v (p,f) = — v x B ■ V v f = —vxB-u = uxB-v 

and by using the fact that C~ l acts only in v, we obtain 

J = j {C^v) ■ (E + u x B) vfidv. 

Since f (C~ 1 v l )v^ fidv = for i ^ j, we can derive the following Ohm's law 

J=-(E + uxB) (3.6) 



where 

i i r 

(£^ 1 v) ■ vfidv. 



1 _ 1 
rj 3 

The derivation of the dynamics of p, u, 9 is very similar to the derivation of the incompressible 
Navier-Stokes- Fourier system (for instance see [1]) except the presence of E, B and g. The 
local conservation laws of mass and momentum of the first equation in ()3.3|) are written as 



e(d t f 1 ,fi)+V x -(f £ 1 ,vfi) = 

e(dtf[,vfi) + Vx^/f.u® vfi) + e 2 {E\ ■ V v (nfi),v) + e{(v x B{) ■ V„(/x 5 f), v) = 
where we have used j(v x B\) ■ V v (fig\)dv = 0. By letting e go to zero, we see that 

V X '(f,v(i) = and V X '(f,v®v[i)=Q 
which yield the divergence-free condition and the Boussinesq relation 

V x -u = and V x (p + 9) = 0. 
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The limiting momentum equation is obtained from 

(d t f!,vfi) + -V x -(f!,v ® vfi) + e{Ef ■ V v (figt),v) + {(v x B{) ■ V v (jisfi),v) = 0. 

(») g w w 

The first two terms, treated in 0], give rise to (a) — > <9 t u and (6) — > u-V x u + V^p — z/Au and 
here the viscosity i/ depends on the collision operator L. See the section 4 in [4] for more 
detail for the treatment of (b). Note that (c) — > 0. And for (d), we see that 

(d) -> J v(v x B)- V v (pg)dv = -J x B. 

Hence, we obtain the following momentum equation 

d t u + u ■ Vu + Vp = i/Au + J x B . (3.7) 

Similarly from the local conservation law of the energy, we can derive the equation for 6. 
Since 

J ^-(v x Bi) ■ V v (p,gi)dv = - j{v x B x ) ■ vpgidv = 0, 

we indeed get the same Fourier equation as in the pure Boltzmann equation 

d t 6 + u • V0 = kA9. (3.8) 

Together with the limiting Maxwell's equations, we deduce that p,u,6, J, E, B satisfy the 
incompressible Navier-Stokes-Fourier-Maxwell system: 

• p + 9 = 0, V-u = 

• dtu + u ■ Vu + Vp = i/Au + J x B 

• d t 9 + u • V6 = kA6 (3.9) 

• E + u x B = r]J 

• no^od t E - V x B = -Ho J, d t B + V x E = 0, V • B = 
This system was studied in |37j . We have established the following result: 

Theorem 3.1. Let [ff , gf, Ef, Bf] be the classical solutions to (|3.3p . Assume that f±, gf, 
Ef, Bf converge strongly to f, g, E, B as e — > 0. Then the limit f has the form f = p + v-u + 

I 1 2 

( 2 — f)^ • Moreover, p,u,6, J, E,B solve the incompressible Navier-Stokes-Fourier-Maxwell 
system (|3.9[) , 
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3.2 Incompressible viscous MHD 

We now discuss the viscous MHD limit. We consider the same scaling as in (|3.2p and in 
addition, we take cq = e so that the electric permittivity is sufficiently small: 

ed t F £ + v ■ V X F £ + (sE E + vx B £ ) ■ V V G £ = \q{F £ , F £ ), 



:i),G- -! r ■ V ,G £ + (— + • V,F e = -_Q(G . h - 



e e 



(3.10) 



p ed t E £ - V x B £ = -no [ v G £ dv, V • B £ = 0, 
d t B £ + V x E £ = 0, V • E £ = - I G £ dv. 

As before, we consider the following perturbation of F £ , G £ , E £ , B £ around \p, 0, 0, 0]: 

F £ = M (l + e/f), G £ = e/ztf, £ £ = eEf , 5 £ = sB\ 

and we suppose that there exist the classical solutions /f, gf, i?f, 5f to (|3.2|) and further 
assume that as e — >• 0, /f — > /, —tg,Ef—} E, B £ — > i?. Following the same procedure as 
in the previous section, we can derive the equations for p,u,9, J, E, B. The only difference 
from the previous system (|3,9p lies in the Maxwell's equations. The last equation in (|3,10p 
forces a = when e — > and the third equation in (|3.10p yields the Amperes law VxB = /ioJ. 
We obtain the following theorem: 

Theorem 3.2. Let F £ = fi(l + eff), G £ = sp,g\, E £ = eEf, B £ = eB\ be the classical 
solutions to (|3.10p . Assume that ff, g\, Ef, Bf converge strongly to f,g,E,B as e — > 0. 

I 1 2 

Then the limit f has the form f = p + v ■ u + (-^j |)# . Moreover, p, u, 6, J, E, B solve the 

following incompressible Navier- Stokes- Fourier-MHD equations: 

• p + 6> = 0, V-u = 

• dtu + u ■ Vu + Vp = uAu + J x B 

• d t 9 + u-V9 = kA6 (3.11) 

• E + u x B = nJ 

• V x B = hqJ, d t B + V x E = 0, V • B = 
3.3 Incompressible inviscid MHD 

To see how the incompressible inviscid MHD system can be derived, we introduce the following 
scaling to (|3.ip by letting t = e 3 t, x = e 2 x, E = -Jj^E, B = -^B, G = -^G. By dropping ~, 

ed t F £ + v ■ V X F £ + {eE £ + vxB £ )- V v G £ = \q(F £ ,F £ ), 

£dtG e + v . VxG e + 1(E £ + VL^L) . VyF e = -Iq( G £ , F £ ), 



p e £d t E £ - 4rV x B £ = / v G £ dv, V • B £ = 0, 

d t B £ + V x E £ = 0, V -E £ = -^r I G £ dv. 

eoe 



(3.12) 
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We consider the following perturbation of F £ ,G £ ,E £ , B £ around [p, 0, 0, 0]: 

F £ = n(l + efl), G £ = epgf, E £ = sEf, B £ = eB{. 
We plug this expansion into (I3.12P to obtain the similar expression to (13 . 3|) : 

ed t fl + v • V x f{ + e(eE{ + v x B\ ) • -V v (pg\) 

= ^-{Q(m/i , A*) + Q(n, l*ff) + eQ(pf!,pf[)} 

ed t gl + v ■ V x g{ + + - * ^ ) ■ -V v (p + e/i/f) 

! £ " (3.13) 
= T~{Q(A*SiiA*) + eQ(M^>^/i)}i 

fi e ed t El -\v xB\ = -^ [ vpg\ dv, V • fl? = 0, 

d t £f + V x £f = 0, V • £f = -L / W f cfo. 

As before, under the convergence assumption of the solutions, we deduce that / = p + v ■ 

u + (•sj |)# . and also we can derive the same Ohm's law (|3.6j) . We can follow the same 

procedure as in the previous sections to derive the equations for p, u, 6, J, E, B. The difference 
is the faster relaxation \ than the previous cases and the diffusion for fluid variables doesn't 
appear at the first order. We can obtain the following theorem. 

Theorem 3.3. Suppose that there exist the classical solutions /f, gf, Ef, Bf to f|3. 13j) and 
further assume that as e — > 0, ff — > f , gf — > g, Ef — > E, B\ — > B. Then f = p + v ■ u + 

I 1 2 

(•sj |)0 . and the limiting equations can be recorded as follows: 

• p + = 0, V • u = 

• dtu + u ■ Vu + Vp = J x B 

• 3 t 9 + u- V0 = O (3.14) 

• E + u x B = rjj 

• V x B = p J, d t B + V x E = 0, V • B = 

which we refer to the incompressible inviscid (Euler) resistive MHD equations. 

We remark that the incompressible, inviscid, ideal MHD equations are valid when the 
resistivity r\ is negligible. 



3.4 Different masses and Hall effect 

We now look at the case when m+ 7^ m_. Moreover, we consider the quasi-neutral regime 
where n + = n_ = 1, we assume the same charges e+ = e_ = 1, and that we choose To = 1 
such that 

J p + dv = J p-dv = 1 
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where fi + and are given in (|1.6p . We also assume that 

(m+/i+ + m-fi-)dv = m+ + m_ = 1. 



We take the same parabolic scaling for t, x, E, B as done in (13.10P to (II. ip except for the 
scaling of G: 

edtF e + „ . VxF e + J_(£* + w x — ) • V„i^ = F%) + -Q(F^, Ff ), 

I B £ 1 1 1 j 

e<9 t Ff + v ■ V x F e _ (E £ + v x — ■ V„Ff = -Q(Fi,F%) + -Q{F £ _,F £ _). 

m,- e e e 

We then consider the following expansion of the solutions near the global Maxwellians /x+ 
and 

i^=// + (l + e /|), Fl + Eft), E £ = eE{, B £ = eB{. (3.16) 

The system (|3. 15|) can be rewritten as 

edtfl + v ■ V x f + + —^—(E{ + vx^)- V„(m + + ^+/+) 



e V/-/ /" 



ed t fi + v • Vs/f —(Ef + vx^-)- V„(m- + e/i_/f 



(3.17) 



where we have used .C 1 and £ 2 to denote each component of the linearized collision operator 
£ given in (jl.7p . By taking e — > 0, we see that 



J 

namely 

f f + (t,x,v)\ I / 1\ _. /0\ (m + v 



m + \v\ 2 _ 3, 
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where n + ,n ,u,8 are to be determined. Then the local conservation laws of mass and 
momentum yield 

•»(©-GD>- +v - <(£)■(>-* 

+ ((eEt + v x BD-V^+fl - t;) = 0. 

As e -)■ 0, since ((v x Bi)-V„(/x+/+ - //_/_),«) = -(v,fi + f + - //_/_) x Bi = x Bi = 0, 
we deduce that 

Vl ' ( (/J'(o) >M = ' Vl ' ( (/J'(o) >M = v,-« = o 

\/-/ \m-V ® v J 
where = n + + n_ . We remark that 

J = y «(//+/+ - fi-f-)dv = 0. (3.18) 

By assuming that 

lim - J £ = Ji (3.19) 

one can derive the MHD type equations containing J\ x Z?i as before. 

To discuss the Ohm's law, we write out the dynamics of the charge density G £ := [i+fl — 
fi-fi from dSUD: 

£dtG e + W . V:E G £ + (E{ + vx^)- V v (^- + ^ + e(^± + ^)) 

e m+ m_ m + m_ 

+ {Q(m+4,m+/+) + Q(m+/+^-/I) - -Q(M-/i,M-/i)} 

Now projecting (|3.20|) onto v, we obtain the dynamics of J e = (vG £ ): 

e9 t J £ iV,-^® vG £ ) - ( + + e(-t- + — 

Av,H+f+) (v,n-fl) 1 (3.21) 



(*) 

To derive the Ohm's law, we split the last term (*) in the left-hand-side of (|3.2ip into two 
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parts: 



(v,m + (i + f^ + m-fj,-fl.) 1 1 
* = ± x B{ + )(v,V+fl - n-fi) x B\ 

1 „ „ , 1 1 x r 

-u e x B\ + ( )J e x Bf . 



(3.22) 



m + m_ m_|_ m_ 



(**) 

By (|3.18p . (**) — > 0. Hence the Ohm's law for J\ = lirn e _>.Q \J e is given by 

r/Ji = — {E x +uxB 1 ) + V X P + C 

in the limit e — >• 0. Here [E\,B{\ = lim. e ^.o[Ef,Bf\, u = lim £ ^.o(v, m + fi + f^_ + m_/i_/i), 

V^-P = — Vjb^w ® vGq) where Go = lim e ^o G e , and the pressure P may depend on ^ 

and (-3- \-)9. Notice that we didn't have the pressure term in the previous simplified 

cases because we started with the strong quasi-neutral scaling where Go = 0. The last term 
C coming from the last term in (|3.2ip is given by (v, Q([j,+f+, //-/-) — Q(fJ,-f-,fJL+f+)). In 
fact, /+ and /_ are purely hydrodynamic described by n^, u, 9 and hence C may depend on 
n^, u, 6 quadratically. Notice that C vanishes when m+ = m_. 

Now we briefly discuss the Hall effect in the above setting. As in the compressible case, 
the Hall effect is not seen in the first order limit. The Hall effect is represented by the Hall 
term J\ x B\, which should be seen from (**) in (|3.22p . Since lim e ^o J e = by (|3.18p . 
(**) is at the leading order. If we approximate J £ by eJ\ by using (|3.19p . we see that 
(**) = e(^p- — ^z)Ji x B\ + o(e). Therefore, we deduce that the Hall term J\ x B\ in (|3.2ip 
can be captured at the next order by including 0(e) terms in the spirit of Chapman-Enskog 
type expansions. It is interesting to point out that (**) is always for m+ = m_ and thus 
it will not affect the Ohm's law at all; for instance, see Section 3.1 and 3.2. This confirms 
that the Hall effect is a two- fluid effect. In the case of m+ ^ m_, the general Ohm's law 
(|3.2ip is much more complicated in that not only the Hall effect but also other higher order 
effects might appear. This motivates us to study the higher order expansion which will be 
investigated in a separate article. 



4 Conclusion 

We would like to conclude the paper with few remarks: 

1) In this paper, we tried to give a formal derivation in the spirit of [1] to MHD type 
models starting form the Vlasov-Maxwell-Boltzmann system. For compressible models, we 
were able to derive the l| Euler-Maxwell system. It seems to us that the classical two fluid 
Euler-Maxwell system can not be derived from the model we started from and we plan to 
come back to this problem in a forthcoming work 

2) We were also able to derive Ohm's law and to compute the resistivity from the Boltz- 
mann kernel. This can be compared to the formulas giving the viscosity v and heat diffusivity 
k in [3]. It is worth pointing out that in our derivation the Hall effect which yields a correc- 
tion to the Ohm's law can only be seen as a higher order correction. This is similar to the 
compressible Navier-Stokes system which is a correction to the compressible Euler system. 

3) There are few other important scaling parameters that we did not take into account in 
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these derivation. One of them is the mass ratio between the electron and the ions. Usually 
it reflects the fact that the electrons attain their equilibrium much faster than the ions. 

4) As observed by Grad [23J, Ohm's law is also valid through an entirely different mecha- 
nism, the so-called gyro-oscillation. We will leave this (magnetic) gyro-effect with or without 
collisions for future study. 

5) After most of this work was done, we learnt that D. Arsenio and L. Saint-Raymond 
[3] are also studying similar limits starting form the renormalized solutions of DiPerna and 
Lions. 
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